Nodeless two-gap superconductivity in stoichiometric iron pnictide LiFeAs 
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The variations of in- and inter- plane London penetration depths, AX(T), were measured using 
a tunnel diode resonator in single crystals of the intrinsic pnictide superconductor LiFeAs. This 
compound appears to be in the clean limit with a residual resistivity of 4 (T — > 0) to 8 (T c ) fiQ-cm 
and RRR of 65 to 35, respectively. The superfluid density, p s (T) = A 2 (0)/A 2 (T), is well described 
by the self-consistent two-gap 7— model. Together with the previous data, our results support the 
universal evolution of the superconducting gap from nodeless to nodal upon departure from optimal 
doping. We also conclude that pairbreaking scattering plays an important role in the deviation of 
the low-temperature behavior of A(T) from exponential in Fe-based compounds. 
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Studies of the superconducting gap structure play an 
important role in the determination of the mechanism re- 
sponsible for superconducting pairing. In iron-based su- 
perconductors [l|, the situation regarding the gap struc- 
ture remains controversial. In the "1111" RFeAs(0,F) 
compounds (R = rare earth) at optimal doping tun- 
neling 0] and angle-resolved photoemission spectroscopy 
(ARPES) 3 experiments have found a full superconduct- 
ing gap. An exponential low-temperature variation of 
the penetration depth, A(T), was reported in compounds 
with magnetic rare earths Smllll [4J and Prllll [5j, but 
the non-magnetic Lallll has shown a close to T 2 power- 
law behavior, incompatible with a clean, full gap Q. A 
similar power-law behavior of A(T) was found in opti- 
mally electron -doped Ba(Fei- x T x ) 2 As2 ("BaT122",T = 
transition metal) lZH10| and hole-doped Bai_ a; K a; Fe2As2 
("BaK122") [H[li]. It was suggested that such a non- 
exponential behavior comes from a full gap with pair- 
breaking scattering 1(| 11, 12- 3- On the other hand, 
in clean (as suggested by observations of quantum oscil- 
lations |19l42l|) isoelectronic P-doped BaFe2(Asi_ a; P a ;)2 
("BaP122") [H, and low-T c stoichiometric LaFePO (H] 
and KFe 2 As 2 ("K122") [3 SBj], the gap appears to be 
nodal, which was suggested to be an intrinsic behav- 
ior of clean Fe-pnictides. However, nodal behavior was 
observed also in dirty systems. Doping-dependent gap 
anisotropy 26, 27] and nodes [jlHU were reported for 122 
crystals. It was suggested that a full gap at optimal dop- 
ing evolves into a 3D nodal structure when doped toward 
the edges of the superconducting "dome". Hereinafter, 
when we refer to a "full gap", we cannot exclude the 
possibility of some angular variation (e.g., see Refs. [In- 
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Since doping inevitably introduces scattering 32 , 
which is pairbreaking in iron pnictides T3, 12, 14, 16- 
HH, measurements of stoichiometric intrinsic supercon- 



361 ] is among very few such compounds. 
It is one of the cleanest systems with a high residual re- 
sistivity ratio (RRR) of about 50 [3611, m uch higher than 
BaP122 (5 to 8 for different doping) [33 , BaCol22 (3 to 
4) [H and BaK122 (7 to 10) pure Bal22 (7 to 10 
under pressure) jipj . though still below pure K122 (over 
1000) 21]. Since T c of LiFeAs decreases with pressure 
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35l |41|, which is observed only in optimally and over- 
doped compounds (HI, we can assign its "equivalent" 
doping level as slightly overdoped, as opposed to under- 
doped NaFeAs, whose T c goes through a maximum with 
pressure [43]] and heavily overdoped K122. This doping 
assignment is consistent with the temperature-dependent 
resistivity, discussed later. With the much reduced effect 
of pairbreaking scattering, comparison of these stoichio- 
metric compounds can bring an insight into the intrinsic 
evolution of the superconducting gap. 

High chemical reactivity leading to quick degradation 
in air makes LiFeAs challenging to study. A single, fully 
isotropic gap 44], as well as two-gap s± pairing [45[ 
were reported from ARPES. Two-gap superconductivity 
is also supported by magnetization [46j, |47j, specific heat 
[IH and nuclear magnetic resonance (NMR) [49| . 

In this Letter, we report on the in-plane London pene- 
tration depth, A(T), in single crystals of LiFeAs measured 
using a tunnel diode resonator (TDR) [Hoj ] . The super- 
fluid density can be well fitted with the self-consistent 
clean two-gap 7— model [5l|. Our results imply that the 
ground state of pnictide superconductors in the clean 
limit at optimal doping is given by s± symmetry with 
two distinct gaps, Ai/T c ~ 2 and A2/T c ~ 1. 

Single crystals of LiFeAs were grown in a sealed tung- 
sten crucible using Bridgcman method sua and were 
transported in sealed ampoules. Immediately after open- 
ing, (0.5 - 1) x (0.5 - 1) x (0.1 - 0.3) mm 3 pieces of the 
same crystal (all surfaces cleaved in Apiezon N grease) 
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TABLE I. Sample properties and parameters of exponential and power-law fits. (See text) 



sample 


r c onset (K) 


T c (K) 


AT C (K) 


ri 


A' (pm/r) 


A (pm/K 3 ' 1 ) 


A {k B T c ) 


#1 


18.2 ±0.1 


17.5 ±0.1 


0.7 ±0.2 


3.39 ± 0.04 


64.8 ±4.8 


106.1 ±0.7 


1.09 ±0.02 


#2 


18.1 ±0.1 


17.2 ± 0.1 


0.9 ±0.2 


2.96 ±0.04 


136.9 ± 10.2 


107.1 ±0.7 


0.94 ±0.02 


#3 


18.0 ±0.1 


16.9 ±0.1 


1.1 ± 0.2 


3.05 ±0.06 


119.7 ± 12.4 


110.1 ±0.9 


0.95 ±0.02 



were used for TDR, transport and magnetization mea- 
surements. Small resistance contacts (~ O.lmfi) were 
tin-soldered [52| and resistivity was measured using a four 
probe technique in Quantum Design PPMS. The transi- 
tion temperature, T c , was determined at the maximum 
of the derivative dAX(T)/dT, TableQ] The London pene- 
tration depth was measured with the TDR technique (for 
review, see (HI). The sample was inserted into a 2 mm 
inner diameter copper coil that produced an rf excitation 
field (at / w 14 MHz) with amplitude H ac ~ 20 mOe, 
much smaller than H c \. Measurements of the in-plane 
penetration depth, AA a h(T), were done with H ac \\ c- 
axis, while with H ac J-c we measured AX c , m i X (T) that 
contains a linear combination of A a & and A c [9]. The 
shift of the resonant frequency, Af(T) = — G4ttx(T), 
where x(T) 18 the differential magnetic susceptibility, 
G = foV s /2V c (l — N) is a constant, N is the demag- 
netization factor, is the sample volume and V c is the 
coil volume. The constant G was determined from the 
full frequency change by physically pulling the sample 
out of the coil. With the characteristic sample size, R, 
Attx = (\/R)ta,nh(R/X) — 1, from which AA can be ob- 
tained [53L [Hi]. 

The main panel in Fig. Q] shows the temperature- 
dependent resistivity, p(T) (left axis), and skin depth, 
S(T) (right axis). p{T) up to room temperature is 
shown in the top inset. The residual resistivity ratio, 
RRR = p(300iT)/p(20 K) = 35 and it reaches the value 
of 65 when extrapolated to T — using a 2 nd order poly- 
nomial. This behavior is consistent with the T-dependent 
resistivity of BaCol22 in the overdoped regime [55j. The 
calculated skin depth, S P (T) — (c/2n)yj p/ f , compares 
well with the TDR data for T > T c , where A/// = 
G[l-n{tanh(aR)/(aR)}], a = {l-i)/5 when we use 
p(300K)—250 pft-cm, the lowest directly measured value 
among our crystals. A very good quantitative match of 
two independent measurements gives us a confidence in 
both resistivity data and the TDR calibration. 

To check for degradation effects, a sample was inten- 
tionally exposed to air for an hour and the measurements 
were repeated, as shown in the lower inset in Fig.[TJ After 
the exposure, the sample surface lost its shiny metallic 
gloss, but the transition temperature and width remained 
the same. At the same time, the total frequency shift 
through the transition (proportional to the sample vol- 
ume) decreased, suggesting that the surface had degraded 
without introducing any change in the bulk. This experi- 
ment clearly supports bulk uniform superconductivity of 
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FIG. 1. Left axis: resistivity (symbols) along with the 2 nd or- 
der polynomial fit from T c to 50 K used to determine the resid- 
ual resistivity, p(0) = 3.7 (i£l-cm. Right axis: skin depth, S, 
measured by TDR, Stdr compared to that calculated from re- 
sistivity, 8 P . Upper inset: p(T) in the full temperature range. 
Lower inset: TDR data for pristine and air-aged samples (see 
text) . 



our samples. 

AA a b(T) in three LiFeAs crystals is shown up to T c /3 
in Fig. [2] by solid dots. AA a &(T) was analyzed using 
(1) power-law, AA(T) = AT 11 (with A and n being free 
parameters) and (2) exponential BCS form, AA(T) = 
Ao"\/7rAo/2Texp (— A /T) (with Ao and A as free pa- 
rameters). The best fit results for sample #1 are shown 
with solid (power-law) and dashed (exponential) lines. 
The fit residuals are shown in the inset. The exponen- 
tial fit quality is as good as the power-law, although 
A /T c = 1.09 ± 0.02 is smaller than the value of 1.76 
expected for a conventional single fully-gapped s-wave 
pairing and Ao = 280 ± 15 nm is somewhat larger than 
the experimental 200 nm 44, Ht]]- This is naturally ex- 
plained by two-gap superconductivity in LiFeAs. The su- 
perconducting T c and best fit parameters (obtained from 
fitting up to T c /3) for all samples are summarized in Ta- 
ble |H T c onsct was defined at 90% of the rf susceptibility 
variation over the transition: the mean T c was defined at 
the maximum of dA\(T)/dT and AT C = T c onsct - T c . 
Ao/T c from the single-gap exponential BCS behavior. 
The power-law coefficient A' was obtained with the ex- 
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FIG. 2. Main panel: AX a b{T) in three LiFeAs crystals (solid 
dots) and A\ c , m i x (T) for sample #2 (crosses). Analysis 
(shown for #1) was done assuming both power-law (solid 
lines) and exponential (dashed line) T— dependences. The 
data for samples #2 and #3 (shifted vertically for clarity by 
10 and 20 nm, respectively) were analyzed in a similar way, 
see Table HI Inset: comparison of the fit residuals for sample 
#1 for the power-law and exponential functions. 



ponent n as a free parameter, while A was obtained with 
a fixed n = 3.1 (average of 3 samples). Crosses in Fig. [5] 
show AX Cj7n i x (T) for sample #2. A clear saturation of 
A\ Ci mi X (T) at low temperatures suggests exponential be- 
havior of A c . 

FigureGUa) shows AA(T) vs. T n with n = 3.1 which is 
the average exponent for three samples. The dependence 
of the parameters, n and A, on the fitting temperature 
range is summarized in Fig. [3] (b) and (c), respectively. 
As expected, the exponent n is more scattered for the 
shortest fit interval, otherwise n and A do not depend 
much on the fitting range from base temperature to 6 K 
and give n > 3 for all samples, with the average value 
3.13 ± 0.23. With n fixed at this average value, we de- 
termined the prefactor A = 107.8 ± 2.1 ^m/K 31 . 

The superfluid density, p s (T) = [1 + AA(T)/A(0)]~ 2 
is the quantity to compare with the calculations for dif- 
ferent gap structures. Figure |H shows Ps(T) for crystal 
#1 calculated with A(0) = 200 nm @,|57j. A noticeable 
positive curvature above T c /2 is similar to other Fe-based 
superconductors [H^] and MgB 2 [H} . This suggests multi- 
gap superconductivity in LiFeAs, which we analyze in 
the framework of the self-consistent 7— model [Hl|. Since 
LiFeAs is a compensated metal and its bandstructure is 
fairly two-dimensional [45j . we may introduce the mass 
ratio on two cylindrical bands, p = mi/m2, whence the 
partial density of states of the first band, n\ — /i/(l + /i). 
The total superfluid density is p s = jpi + (1 — 7)^2 with 
7 = 1 / (1 + We also use the Debye temperature of 240 
K [48j to calculate T c , which allows fixing one of the in- 
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FIG. 3. (a) AA vs. T 3 ' 1 (average exponent n over 3 samples); 
(b) exponent n and (c) prefactor A, obtained by fitting to 
AA(T) oc AT n for various upper temperature limits shown 
on the :r-axis. 



teraction constants, An. This leaves three free fit param- 
eters: the in-band, A22, and interband, A12, interaction 
and the mass ratio, /i. Figure 0] shows that p s (T) can 
be well described in the entire temperature range by this 
clean-limit weak-coupling BCS model. In the fitting, the 
two gaps were calculated self-consistently (which is the 
major difference between this one and the popular, but 
not self-consistent, a - model 60]) and the self-consistent 
Ai(T) and A2(T) are shown in the upper inset in Fig. 21 
while the gap ratio is shown in the lower inset indicat- 
ing strong non-single-gap-BCS behavior of the small gap. 
The best fit, using Levenberg-Marquardt nonlinear sim- 
plex optimization in Matlab, gives Ai(0)/T c ~ 1.885 and 
A2(0)/T c ~ 1.111. As expected, one of the gaps is larger 
and the other is smaller than the single-gap value of 
1.76, which is always the case for a self-consistent two- 
gap solution. The best fit parameters are: An = 0.630, 
A 22 = 0.642, A12 = 0.061 and \i = 1.384. The deter- 
mined mass ratio gives n\ = 0.581 and 7 = 0.419. This 
is consistent with bandstructure calculations that yield 
rii = 0.57 and p = 1.34 6l|, and ARPES experiments 
that find fi m 1.7 [iij . The effective coupling stren gth , 
X e ff = 0.374, is not far from 0.35 estimated for 122 [62j 
and 0.21 for 1111 [63} pnictides. The electron band with 
a smaller gap gives about 1.5 times larger contribution to 
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FIG. 4. Symbols: superfluid density, ps{T) calculated with 
A(0) = 200 nm. Solid lines represent the fit to a two-gap 
7— model, p s = 7P1 + (1 — 7)p2- Dashed line is a single-gap 
BCS solution. Upper inset: superconducting gaps, Ai(T)and 
A2(T) calculated self-consistently during the fitting. Lower 
inset: A1/A2 as a function of temperature. 



the total p s resulting in a crossing of the partial densities 
at low temperatures. Similar result was obtained from 

We stress that while 



46|. 



magnetization measurements 
Ai(T), A2(T) and p (hence, rii, n-i and 7) and A e // are 
unique self-consistent solutions describing the data, the 
coupling matrix \j is not unique. There are other com- 
binations that could produce similar results, so it seems 
that Xij has to be calculated from first principles [6l| . 

In conclusion, we find that in the clean limit, 
optimally-doped Fe-based superconductors are fully 
gapped, but most measurements are affected by pair- 
breaking scattering 

0,1 

Til iRj - fljj . This conclusion is in 
line with studies of thermal conductivity (which is not so 
sensitive to scattering) in BaCol22 [271 12811 . On the other 
hand, intrinsic K122 reveals a nodal gap I24j.l25l . which is 
also found in heavily overdopcd BaCol22 [28]. Overall, 
this establishes a common trend for all Fe-based super- 
conductors to have a superconducting gap that evolves 
from full to nodal when moving toward the edge of the 
superconducting dome. 
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